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We show that the Bethe-Peierls (BP) boundary condition should be modified for ultracold atoms 
with spin-orbit (SO) coupling. Moreover, we derive a general form of the modified BP boundary 
condition, which is applicable to a system with arbitrary kind of SO coupling. In the modified 
BP condition, an anisotropic term appears and the inter-atomic scattering length is normally SO- 
coupling dependent. For the special system in the current experiments, however, it can be proved 
that the scattering length is SO-coupling independent, and takes the same value as in the case 
without SO coupling. Our result is helpful for the study of both few-body and many-body physics 
in SO-coupled ultracold gases. 
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I. INTRODUCTION 



In the study of ultracold atomic gases the 
Bethe-Peierls (BP) boundary condition [l[ is widely used 
as a replacement of realistic interaction between two 
atoms. With this approach, one only needs to solve the 
Schrodingcr equation with the Hamiltonian for atomic 
free motion, and thus the calculation is significantly sim- 
plified. As a result, the BP boundary condition is very 
useful in the research of few-body and many-body physics 
in ultracold gases, especially those with large interatomic 
scattering lengths. Many achievements have been ob- 
tained. For instance, for two-component Fermi gases, 
Petrov et. al. obtained the atom-dimer [f| and the 
dimer-dimer [(| scattering lengths, and Werner et.al. Q 
redcrived the well knonwn Tan's relations using 
the BP boundary condition. 

In recent years, a class of synthetic gauge fields and 
spin-orbit (SO) c oupling has been realized in ultracold 
Bose gases [14l - [20j and degenerate Fermi gases [2l|, [22| ■ A 
considerable amount of theoretical interest has been stim- 
ulated to understand the SO-coupling effect in both few- 
body [23rl27j and many-body physics [H-HH, inclu ding 
the gases with large interatomic scattering lengths |35r - 
H3|. It becomes now an urgent task to carefully exam- 
ine the BP boundary condition in SO-coupled ultracold 
gases. 



II. MODIFIED BP BOUNDARY CONDITION 
FOR SPIN-I/2 FERMIONIC ATOMS 

In this section we shall consider a system of two 
spin-1/2 fermonic atoms with a short-ranged and spin- 
dependent interaction potential U(r), where f= (x,y,z) 
is the relative position of the two atoms. The interaction 
U(r) has an effective range such that U (f) ~ for 
r = \r\ > ?v Furthermore, we shall focus on the case of 
low-energy scattering for which the difference e between 
the energy of atomic relative motion and the scattering 
threshold is much smaller than 1/r 2 . 

Let |t) and \\.) represent the spin eigen-states of a single 
atom, the quantum state of the relative atomic motion 
can be described by a spinor wave function |i/>(r)): 



|^(r)> = TPs (r) |S) 



3=1 



IT, 



(1) 



where |S) = (It) ill) 2 - ||)i|t>2)A/2 is the singlet spin 
state and |Tj) (j = 1, 2, 3) are the three triplet states. In 
dilute ultracold gases, the interatomic distance is much 
larger than the effective range r* , and the physical prop- 
erty of the system is determined by the behavior of the 
wave function \tp(r)) in the region r > r*. 

Our task is to investigate the behavior of the wave 
function \ip(f)} in the presence of SO coupling, and then 
establish the correct BP boundary condition. 



In this paper we show that in SO-coupled systems, the 
BP boundary condition should be modified and moreover 
derive a general form of the modified BP boundary con- 
dition that is applicable to a system with any kind of SO 
coupling and arbitrary atomic spin. The relevance to the 
current experiments is discussed. 



A. Without SO coupling 

For completeness, we start with the case without SO 
coupling, for which the relative motion of the two atoms 
is governed by the Hamiltonian 

H=p 2 + B + U(r) , (2) 
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where p = — tV is the relative momentum and the nat- 
ural units h = m = I (m is the single-atom mass) are 
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used. Operator B acts in the spin space and accounts for 
the possible r- independent contribution, e.g., from the 
Zeeman effect. We assume the difference of the eigen- 
energies of B are much smaller than 1 jr 1 . 

We first consider the property of \^{r)) in the short- 
range region < r << \Jyfe. According to the low- 
energy scattering theory (appendix A), when \ip(r)) is a 
low-energy eigenf unction of H , one has 

|V(r)) oc Q- - ^ |S) for r* <r « 1/Ve, (3) 

with the scattering length a being determined by the de- 
tail of U(r). Note that Eq. ([3]) does not depend on the 
eigen- value of H for \ip(r)), and is thus applicable to to 
all low-energy wave functions. 

With Eq. ([3]), one can obtain the behavior of a low- 



> 



r*. 



energy wave function \ip(r)) in the whole region r 
Let \4>{r)) be the solution of the Schrodinger equation 
with Hamiltonian p 2 together with the BP boundary 
condition 



lim \4>{r)) oc 
>o \ r 



1 1 



|S)+0(r) 



(4) 



the realistic wave function \tft(f)} and the pseudo one 
|c/)(r)) will have the same behavior for r > r*. There- 
fore, in the study of few-body and many-body physics in 
systems with interaction U{r), one can replace the real- 
istic potential U(r) by the BP condition Theoretical 
calculations can be greatly simplified. 



where B{K) = + + A(<ri 1) + &P)K x /2 and the 
c-numbcr K = (K x , K y , K z ) is just the total momentum 
of the two atoms. 

For the aim of establishing a correct BP boundary con- 
dition, one should also examine the behavior of the eigen- 
function \tp(r)) of H in Eq. ([5]) in the short-range region, 
now defined as r, < r << r s with r s = min(l/v / £~, 1/A). 
Comparing Eq. ^ to Eq. ©, one finds that due to the 
SO coupling, the Hamiltonian H is modified by a term 

A(a'i 1 ^ — o)f^)Vx- This term exists in the whole range of 
the interatomic distance r, including the short-range re- 
gion and the region r < r». Therefore, the short-range 
behavior of \tp(r*j) can no longer be described by Eq. (J3J, 
and the BP boundary condition in Eq. §4§ cannot be di- 
rectly applied. 

To overcome this difficulty, we introduce a unitary 
transformation (rotation) lZ(r) as 



1l(r) = e iX & l) -*? ) >/\ 



(7) 



with x the relative position in the x direction, and define 
the rotated wave function \ip(r)}-n as 



\^{r)) K = 1l{r)\iP{r)). 



(8) 



An immediate observation is that since \ift(f)) is an eigen- 
function of Hamiltonian H, the rotated wave function 
\ 1 P(f))R is an eigenf unction of the rotated Hamiltonian 



H R = K{r)HK\r). 
Straightforward calculations yield 



(9) 



B. With one-dimensional SO coupling 



H R =p 2 + W(f) + U R (f) 



(10) 



We now consider a simple case with one-dimensional 
SO coupling. Without loss of generality, the single-atom 
Hamiltonian of the system can be written as 



p2 

Hib = — + Xcr z P x 



Z, 



(5) 



where P is the atomic momentum, a is the Pauli op- 
erator, A indicates the intensity of the SO coupling, 
and Z accounts for the residual spin-dependent part. 
This paper shall consider the case of weak SO coupling 
A << 1/r*, as in the current experiments. We assume the 
difference between the eigenvalues of Z is much smaller 
than l/r|. 

The total Hamiltonian of the two atoms is given by 



( 2) 

H lh ' + U(f), where is for the ith atom 



(0 



Because the total momentum of the two atoms is con- 
served, the relative motion can be separated from the 
mass-center motion. The Hamiltonian for the relative 
motion is then 



H = p z +X (a 



Px + B{K) + U{r) = H a + U(r), 

(6) 



with 



U R (f) = K(f)U(r)K*(f). 



(11) 



W(r) = n(r)B(K)nHr)--(a^-a^r.(l2) 

Equation (|10p shows that the SO coupling disappears 
in the rotated Hamiltonian H R . Furthermore, in the re- 
gion r « 1/A, we have W(r) ~ W(0), and then 



SR 



f +W(0) + U K (r). 



(13) 



Thus, the eigenfunctions of H R and -ffsR have the same 
behavior in the short-range region < r << r s . Fur- 
ther, Eq. {jT3j) has the same as Eq. ([2]). Therefore, the 
eigenfunction \ip(r))n behaves as 

h/>(r)>R oc f - - —) |S) for r,<r«r s , (14) 
\r a R J 

analogous to Eq. ([3]). It should be pointed out that a R 
is the scattering length with respect to the rotated in- 
teraction potential Un(r). On the basis of Eq. ([H)l . the 
behavior of the wave function {^{f 4 )) in the unrotatcd 
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frame can be obtained by the inverse unitary transfor- 
mation It) (r) . The result is 

\m) 4- -) is> - i~{t z i - o z2 ) ■ (-) is) , 

\r a R J 2 \r/ 

for r * ^5 r « r s- (15) 

As in the above subsection, Eq. (ITSl does not depend on 
whether or not \ip(r)) is an eigenfuction of H , and thus is 
generally applicable. Finally, we emphasis that the last 
term in Eq. (| X5[) is in the order of unit and cannot be 
neglected. 

Let \4>{r)) be the wave function given by the 
Schrodingcr equation with Hamiltonian Hq in Eq. 
together with the modified BP boundary condition 

lim 10(F)) cc ( - - -L ) ISM^W-ai 2 )).^) |S>+0(r). 
r->o \r a R J I \r 1 

^ (16) 

It is clear that in the whole region r > r*, \4>(r)) has 
the same behavior with the solution of the Schrodingcr 
equation with Hamiltonian H in Eq. (j6j). Therefore, the- 
oretical calculation can be done by replacing U (F) with 
condition (fT6"|) . 

So far we have obtained the the modified BP boundary 
condition (|16[) for the system with one-dimensional SO 
coupling. Comparing Eq. (fT6|) to Eq. (j4]), we find that the 
SO coupling has two effects on the BP condition. First, 
the modified BP condition includes an anisotropic term 
— iX(a z ^ — a^)x/(2r)\S). Second, the scattering length 
a R is determined by the detail of the rotated interaction 
potential C/r.(F) defined in Eq. (fTT|) . 

Relevance to the current experiments. In general, to 
obtain the value of a R we need to explicitly solve the 
Schrodingcr equation with potential £/r(F) for a given SO 
coupling. As shown below, however, in the special case 
of the current experiments [U [22| , the SO coupling does 
not change the scattering length. Namely, a R equals to 
the scattering length ao for the case without SO coupling. 

In Refs. [2l], [22(, the spin states | f) an d | I) are two 
hypcrfine states of 6 Li or 40 K atom, and the single-atom 
Hamiltonian in the Schrodingcr picture is given by 

P 2 O r) 

r Si , 2ik r X , ~ -2ik r X\ , u ~ / 17 \ 
#ibs = -y + 2 y (J + e ) + 2 (Tz ' ( 17 > 

with X the single-atom coordinate in the x direction, 
<7_l_ = | fK-i- | and <7_ = a\. Here £1 and 8 are the 
Rabi frequency and the two-photon detuning, respec- 
tively. Then, in the Schrodingcr picture the two-atom 

Hamiltonian is i?2bs = H\hs + -^lbS + ^o(^j where 
the bare inter-atomic interaction potential Uo(r) has the 
scattering length ao- 

The SO coupling term emerges after the spin ro- 
tation along the z-axis or the unitary transformation 
T(X) = exp [—ik r Xa z ] is applied. In the transformed 
picture, the single-atom Hamiltonian is given by Hyo = 
T(X)HibsT^ (X). It can be shown that H\\, takes the 
form in Eq. (J5J) with A = 2k r and Z = Qa x /2 + Sa z /2. 



The two-atom Hamiltonian in the transformed picture 

can be written as H-j^ +h[^ +U(r), and the transformed 
interaction potential U(r) is given by 

U(r) = T(X 1 )T{X 2 )U {r)T\X 2 )T^(X 2 ) 

= e - lk ^ ) ~ a ^^ x Uo{r)e lk ^ ) - a( "^ x , (18) 

where the relative coordinate x satisfies x = X\ — X 2 
and we have used the fact [£/o(F), <7i +o'i 2 ' ) ] = 0, arising 
from the conservation of the total z-component of hyper- 
fine spin during the collision process. We emphasis that, 
in the transformed picture where the SO-coupling term 
a z P x appears, the inter-atomic interaction is not the bare 
potcntianl C/o(F), but the transformed one U(r). 

The scattering length a R in the modified BP bound- 
ary condition (|16[) is determined by Un(r) in Eq. pip . 
Substituting Eq. (fT5|) into Eq. (fTT|) . we find that the ro- 
tated potential J7r(F) reduces to the bare potential-i.e., 
Un(r) = Uo(r). Therefore, the scattering length remains 
unchanged with SO coupling-i.e., a R = ao. 

C. With arbitrary type of SO coupling 

We shall now extend the above treatment to the spin- 
1/2 fermonic system with arbitrary type of SO coupling. 
In this case, the single-atom Hamiltonian can be gener- 
ally written as 

P 2 -, 
H lh = — +XM-P + Z, (19) 

with M and Z are operators in the spin space and the 
maximum eigenvalue of M is of the order of unit. After 
being separated from the mass-center motion, the the rel- 
ative motion of the two atoms is described by the Hamil- 
tonian 

H = p 2 + Xc-p + B(K) = H Q + U(r). (20) 

with c = MW-MP) andB(A') = Z^ + Z^ + X(M^ + 
MW)-K/2. 

As in the above subsection, to investigate the short- 
range behavior of the eigenfunction \ijj(r)) of H , we in- 
troduce a unitary transformation lZ(f) as 

7£(F) = e i^c as x/2 e i\c y y/2 e i\c z z/2^ ^1) 

with c = (c Xl Cy, c z ). The rotated Hamiltonian H R = 
1Z(r)H'R) (r) can be calculated as 

H R =p 2 - 2Xd(Xf) ■ p + W(r) + U R (r) (22) 

with operators d = (d x ,d y ,d z ) and W given by 

d x {Xr) = 0; (23) 
d v (Xf) = e <Ac,»/2^f e -»Ac,«/a _ft(^ftt(r). (24) 

d z (Xr) = C ^-K(r)^KHr), (25) 
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and 



W{r) 



iX 



X' 



W-d(Xf) +K(f)B(K)K^(r) + (26) 
\d{Xr)¥-K{r)^K\r)\. (27) 



Here wc have [/r(t) = lZ(f)U (r)ltf (f) as before. 

Unlike Eq. (fTO)) in the above subsection, the SO 
coupling still exists in the rotated Hamiltonian Hn in 
Eq. ([22]). Nevertheless, according to Eqs. (|23ti25|) . we 
have d(Xf) = 0(Xr). Namely, its zeroth-order contribu- 
tion in i?R vanishes. This leads to the following impor- 
tant property of the eigenfunction |?/>(r))R of H^: 

|V(f)>R cx (- - — ) |S) (for r, < r « r s ) . (28) 



The result in Eq. ([28)) is proved as follows. In the region 
r > r* where Un(r) is negligible, the eigen-equation of 
reads 



p 42 - 2Ad(Ar) • p + IU(r) |^(r)) R = E\i>(r)) R . (29) 



As shown in appendix B, in this region the wave function 
|'0(r))R can be expressed as 

W»0>R = -^|S> + £C n r n |S) 



n=0 



E E £r^,m,(^0)l^ 



I 



(30) 



1=1 mi= — l n=0 



in the spherical coordinate (r,8,<fi). Here Yi mi (8,(f>) are 
the spherical harmonic functions, C„ (n = —1,0, 1, ...) is 
the coefficient of term r n |S), and |^4/, mi , n ) (n = 0, 1, ...) is 
the spin-state with respect to r n Y^ mi (9, (f). Substituting 
Eq. (|3T)|) into Eq. (|2"!)1) and comparing the coefficient of 
the term r~ 2 in both sides, we find that because d(Xr) = 
0(Ar), one has \Ai^ mi ^) = 0. Therefore, in the short- 
range region \tp(r)) R behaves as in Eq. (f2"5|). with the 
scattering length a R determined by both the potential 
Us.{f) and the operator d(Xf). 

With Eq. (|28|) and following the procedure in the above 
subsection, we obtain the modified BP boundary condi- 
tion for the general type of SO coupling as 

Mm Wr)) oc(l--) \S)-iU (-) \S) + 0(r). (31) 
r->0 \r a R J 2 \r J 

As in the above subsection, the value of the scattering 
length a R in general depends on the SO coupling. This 
dependence is also shown in Fig. 3 of Ref. [23( with a 
simple model where the potential U(f) is modeled as a 
spin-independent spherical square well. 



III. MODIFIED BP BOUNDARY CONDITION 
FOR ATOMS WITH ARBITRARY SPIN 



and arbitrary spin. The Hamiltonian for the single-atom 
motion and the relative motion of the two atoms are still 
given by Eqs. (fP9"|) and (f2TJ)) . respectively. 

For simplicity, we first consider the case that the inter- 
atomic interaction U (f) (with scattering length a) is in- 
dependent of the atomic spin. In this case, it can be 
shown that without SO coupling, the low-energy eigen- 
state \i(>(r)) of the relative-motion Hamiltonian (—V 2 + 
U{r) + ZW + Z^) behaves as 



[;--) 



for r* < r « 



(32) 



This is very similar to Eq. ([2]) , but now the r-independent 
spin state \x) is not unique. Instead, \x) can be different 
for different cigen-states \ip(r)). 

In the presence of SO coupling, the short-range behav- 
ior of the eigenfunction \tp(r)) can be obtained via the 
same approach using the unitary transformation lZ(r) in 
Eq. ([2~Tj) . In particular, in the region r << r S) it is suf- 
ficient to keep the lowest-order terms of d(Xf) and W(r) 
defined in Eqs. (|23H26|) . Thus, the rotated wave function 
l^(f))R = ^{^\4>(f)) satisfies the equation 

[f - 2Xg(Xr) -p+ W(0) + U{r)} |^(r)) R = E\^{r)) K , 

(r < < r.) (33) 

where E is the eigen-energy. Here we have used d(0) = 
and the fact Un(r) = U(r) which is because U is spin- 
independent. The operator g = (g x , g y , g z ) is defined as 
ffi(Ar) = r ■ [V di(Xr)\f=o] with i = x,y,z. In Eq. lf33]) 
the term —2Xg ■ p couples the s-wave and d-wave compo- 
nents of \ijj(r))-n. In the Schorcdingcr equation, the cou- 
pling terms are either independent of r or proportional to 
r(d/dr), and thus do not decrease the power of r in the 
wave function \ip(r))n. The estimation with the semi- 
classical approximation d\ip(f))-fi/ dr < \J — U(r)\ip{f)) R 
shows that, in the short-range region r* < r << r s , the 
intensity of this coupling is much smaller than the cen- 
trifugal potential 6/r 2 , which is the energy gap between 
the s-wave and d-wave channels. In addition, for many 
systems this intensity is also much smaller than 6/r 2 even 
when r < r* . An example is a system with Lennard- 
Jones potential U(r) = —c^/r 6 + c^jr 12 . Therefore, for 
these systems we can neglect the SO coupling in the en- 
tire region r << r s . Then one has 

h/>(7)>R oc (^ - ^ |x) (for r* < r « r s ). (34) 

Note that a is still the scattering length of the potential 
U(r). Accordingly, we have the modified BP boundary 
condition 

hm 10(f)) oc (- - -) \ X )-i^c- (?) \ X ) + 0{r). (35) 



Finally, we consider the general case: a system of two 
fermonic or bosonic atoms with any kind of SO coupling 



The situation becomes more sophisticated if U(r) is 
spin-dependent or the SO coupling cannot be neglected 
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in Hr when r < r*. In these cases, 1/a in the modified 
BP boundary condition ([33)) should be replaced by an 
operator in the spin space, which is also determined 
by U-R,(r) and d(\r). The detail is given in appendix C. 

We conclude this section by poin ting out that as in 
Sec. II, in the current experiments fl4T - r20| for bosonic 
atoms with one-dimensional SO coupling, the rotated po- 
tential J7r is equivalent to the bare potential Uq in the 
Schrodinger picture, and the operator d is zero. Then 
the operator in the modified BP boundary condi- 
tion is independent of the SO coupling. For instance, 
for the ultracold gases with spin-1 87 Rb atoms, we have 
An = l/a,QpF=o + l/a2'£ , F=2, where ao (0,2) is the scat- 
tering length with respect to the total atomic spin F = 
(F = 2) and Vf=o,2 are the relevant projection operators. 

IV. DISCUSSION 

In this paper we derive the modified BP boundary con- 
dition for ultracold atomic gases with SO coupling. It 
is shown that the SO coupling brings a new anisotropic 
term to the BP boundary condition, and may change the 
value of atomic scattering length. 

Our result can be used for the research of both few- 
body and many-body problems in SO-coupled ultracold 
gases. For instance, for N spin-1/2 fermonic atoms with 
the Hamiltonian 

N N N 

i— 1 i— 1 i<j 

N 

= HF + ^Ufa), (36) 

i<j 

one can replace the interaction potential U (fij ) by the 
modified BP boundary condition 

lim (ry|$) 
\rij\->0 

+0(r y ). (37) 

In Eq. ([55)1 is the trap potential for the zth atom 

and fij is the relative position between the «th and jth 
atoms; in Eq. ([37)) |4>) is the TV-atom state, \fij) is the 
cigen-state of the relative motion of the (i, j) pair, |S)y 
is the singlet spin state for the two atoms and |$') is a 
quantum state for other atoms. The limit in Eq. (|37[) is 
taken for fixing the positions of other atoms as well as 
the mass center of the (i,j) pair. In the region > r*, 
the solution of the Schrodinger equation with the free 
Hamiltonian Hp under the boundary condition ([3"7| has 
the same behavior as the soulution of the Schrodinger 
equation with the total Hamiltonian Ht- 

Our result is also useful in the current experiments of 
ultracold gases with one-dimensional SO coupling and far 



away from the Feshbach resonance point. As shown in 
Sees. II and III, the scattering lengths in these systems 
are independent of the SO coupling, and thus all the 
terms in the modified BP boundary condition can be fully 
determined with the known parameters. 
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Appendix A: Short-range behavior of scattering 
wave function 

In this appendix we prove Eq. §5§ for the short-range 
behavior of the wave function \ip(r)) in the cases without 
SO coupling. Without loss of generality, here we consider 
the case that \ip(r)) is the scattering wave function and 
then satisfies the Lippmman-Schwinger equation 

|V>(f)) = |^ (0) (^))+ / dr'goiE^r'Mr'me')}, (Al) 

Here E is the eigen-energy of H in Eq. @ with re- 
spect to \i()(r\), \il>(°'{r)) is the incident state and satisfies 
(p 2 + B)|V' (0 '(r)) = E\ipW (r)) . For our system with two 
fermonic atoms, \^°'(f)) is anti-symnetric with respect 
to the permutation of the two atoms. In Eq. (jAl[) the 
Green's operator go(E,r,r l ) is defined as 

e Wn~e n \r-r'\ 

= -£ p-^i—|n><»l (A2) 

n 1 1 

with e n and |n) the nth eigen- value and eigen-state of 
the operator B, respectively. 

Since the potential U(r) is negligible in the region r > 
r*, the integration in Eq. (jAip is only effective in the 
region r' < r*. In the low-energy cases, when r — > 00 
and r' < r*, the function go{E, f, f') becomes very steady 
with respect to r' and we have go{E, r, r ') m go(E, r, 0). 
Therefore, in the limit r — >• 00, the solution of Eq. (|A1|) 
takes the form 

|^(r)) = \^°\r))+g (E,r,0)\x), (A3) 

where the spin state \x) is related to \ip{r)) via the equa- 
tion \x) — J dr'U(r')\ip(r')). Furthermore, due to the 
facts Pi 2 \^(f)) = -\ip(r)) and Pi 2 U(r)P l2 = U{f) with 
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Pi 2 the permutation operator of the two atoms, one finds 
that P 12 U(r)\i/j(r)) = -U(r)\if)(f)). This result yields 



the incident wave function |i/;( )(f)), which is an eigen- 
function of Hq defined in Eq. (f^Uj) . takes the form 



[X> = |S> / df'U(r') (S^(f')}- 



(A4) 



On the other hand, since \tp(r)) is an eigen-state of H 
and the potential U (f) is negligible in the region r > , 
in such a region the wave function \ip{r)} satisfies the 
equation 



|V (0) (r)> = 



2(2tt) 3 / : 



■(l-P 12 y k - r \a(k)} 



(B3) 



with P\2 the permutation operator for both the spin and 
the spatial motion of the two atoms. Eq. (|B3[) leads to 
the result that 



|V- (0) (r)> =O(r ). 



(B4) 



[P 



B)\Hr))=E\4>{?)). 



(A5) 



Therefore, the behavior of the wave function \tp(r)} in the 
region r > r» is determined by Eq. (|A5[) and the bound- 
ary condition (|A3I) in the limit r — > oo. Considering Eq. 
l4| . one can easily prove that the function 



Now we consider the expansion of the Green's function 
g(E,r,Q). Using the fact 



5{r-r') 



ik ^ " ! } ( £ Hk))(a(k)\] , (B5) 



\tp(f}} = \^°\f)) +A a g Q (E,f,0)\S) (A6) it is easy to show that 



with Ao a constant satisfies both of the two conditions. 
Therefore, \ip(r)) satisfies Eq. (|A6[) in the whole region 
of r > . 

To obtain the short-range behavior of \^(r)), one can 
expand Eq. (|A6[) as a series of r, and then neglect the 
high-order terms. Using Eq. (|A2[) and the fact that 
Pi2\'4 )t ^ > \'r)) = — we immediately gets the re- 

sult in Eq. ©: 

Mr)) oc f - - - ) IS) forr*<r<<l/Vi. (A7) 
V r a / 



Appendix B: Proof of Eq. ([30)) 

Now we prove Eq. (f30| for the behavior of |i/>(r))R 
in the region r > r*. To this end, we first consider the 
behavior of the un-rotated eigenfunction IV'(r')) of Hamil- 
tonian H defined in Eq. (|2T)|) . Without loss of generality, 
here we consider the case that \ip(r)} is the scattering 
wave function. Using the approach in appendix A, we 
can prove that when r > we have 



\TP(r)) = \^(7))+A o g(E,r,0)\S) 



(Bl) 



Here E is the eigen-energy of H with respect to \ip(r)}, 
\tp^°\r)) is the incident state and satisfies H \tp^(r)) = 
E\ip(°l(r)) with H defined in Eq. (20]). The Green's 
operator g(E,r,r') is defined as 



g(r],f,r') = 



■8(f-f'). 



(B2) 



Now we expand the r.h.s. of Eq. (|B1|) as a power 
series of r. To this end, we first consider the operator 
F(k) = Ac • k + B(K), with k a constant operator and 
A, c and B(K) defined in Sec. II. For each given vector k, 
F(k) is an operator in the 4-dimensional spin space. We 
denote the a-th (a = 1,2,3,4) eigen-energy and eigen- 
state F(k) as £ (a, k) and \a(k)), respectively. Therefore, 



g{E,f,0) = Y^ / dk 



■ik-r 



\a(k))(a(k)\ 



(2ttY E + i0 + 



k 2 +£{a,k) 



(B6) 

Eq. (|B6j) and the completeness relationship 

J2 a l a (^))( a (^)l = 1 l ea, d to the result 



ik-r 



1 



g(E,r,0) = / dk 

(2tt) 3 E + i0+ - k 2 



e ik-r 
1 



|a(fc))(a(fc)| x 



1 



E + i0+ - 



k 2 + £(a,k) 



E + iQ^ 



(B7) 



It is pointed out that, in the limit r — > 0, the integration 
in the r.h.s of Eq. (|B7[) converges to a constant operator 
in the spin space. On the other hand, we also have 



dk- 



1 



iVEr 



(2tt) j E + iO+ - k 2 



(B8) 



Due to these facts, we have g(E,r,0) oc 1/r + O(r ). 
Substituting this result and Eq. (|B4"|) into Eq. ([BT]) 
and using the relation | - 0(O)r = ^•( 7 ^I ? A( 7? )) with lZ(r) 
defined in Eq. (f2"Tj) , we can find that 



1 



^(r)) R cx-|S) + O(r ), 
r 

Namely, \ip(f))-[i takes the form of Eq. (f3"U)) . 



(B9) 



Appendix C: The modified BP boundary condition 
for atoms with spin-dependent interaction 

In Sec. Ill of our maintext, we derive the modified BP 
boundary condition for atoms with arbitrary spin and SO 



7 



coupling. Our result in Eq. (|35j) is based on the following 
two assumptions: (a) the inter-atomic interaction is spin- 
independent; (b) in the rotated frame, the influence of 
the SO coupling or the term —2g ■ p is negligible in the 
region r < . In this appendix we go beyond these two 
assumptions and derive the general type of the modified 
BP boundary condition for atoms with with arbitrary 
spin and SO coupling. 

We first go beyond the assumption (a) and consider 
the case of atoms with spin-dependent interaction U (r) . 
When there is no SO coupling, the eigenfunction \ijj(r)) 
of the two-atom relative Hamiltonian satisfies 



Z 



(i) 



A2) 



\^{r))=E\^{r)). (CI) 



We assume the spin space of the two atoms is n dimen- 
sional, and denote the eigen-states of Z^> + as \j) 
(j = 1, n). We further define 

- r#(r)> • 

Then IvE'(r)) satisfies the boundary condition |^(0)) = 0. 

We first consider that U{r) is spherical. Thus, ^(r)) 
can be written as 



TO)=5>;(r)|j) 



(C2) 



Then, Eq. (jCip can be re-expressed as the equation for 
We define |$( Q )(r)) as the s-wave solution of this 



equation, with component (r) satisfying the bound 
ary conditions (0) = and 



dr 



$(«)( r ) 



1, for a = j 
0, for a =/= j 



(C3) 



Therefore, the states \(j> {a) (r)) = \& a \r))/r are n special 
solutions of Eq. (|C2[) . In the short-range region, the low- 
energy wave function \4> (f)} behaves as 



\<t>W(f)) = hM a ) 



\T a ) 



(C4) 



with \M a ) and \T a ) are states of atomic spin. 

Furthermore, any s-wave solution \ip(r)) of Eq. (|C1[) 
can be written as the linear combination of (r)), and 
then expressed as 



b a 



-\M a ) 
r 



\T a ) 



for r* < r << r s 



(C5) 

In addition, the low-energy solutions of Eq. (|C1[) with 
high partial waves are negligible in the short-range re- 
gion. Therefore, Eq. (|CJ5[) is actually satisfied by all the 
low-energy solutions of Eq. (jCll) . 

When the states \M a ) with different a are linearly in- 
dependent of each other, we can define an operator A 
which satisfies A\M a ) = \T a ) (in particular, when the 
interaction U is independent of the atomic spin, we have 
A = 1/a). With this definition, the behavior (|C5|) of 
\ijj(r)} can be re- written as 



[--A] 



As in Sec. Ill, the r-independent state |x) in the spin 
space is not unique. Finally, it can be proved that in 
the low-energy limit the above result is also correct when 
U(r) becomes anisotropic. 

In the presence of SO coupling, with Eq. (|C6|) and the 
approach in our maintext we can obtain the modified BP 
boundary condition: 



lim |<Kr)) ex (I - Ab\ \ X )-i^c- ( - ) 



with the operator Ar determined by both the potential 
U^i(r) and the operator d(Xf). 

Finally, if we go beyond the assumption (b) and con- 
sider the case that the SO coupling cannot be neglected 
when r < , we can also follow the above approach, and 
obtain the modified BP boundary condition which has 
the form in Eq. (fCTf . 
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